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We examine the competition and relationship between an antiferromagnetic (AF) 
Mott insulating state and a d x 2_ y 2 superconducting (SC) state in two dimensions 
using semi-analytical, i. e. diagrammatic calculations of the t—U—W model. The AF 
Mott insulator is described by the ground state of the half-filled Hubbard model on a 
square lattice with on-site Coulomb repulsion U and nearest neighbor single-particle 
hopping t. To this model, an extra term W is added, which depends upon the square 
of the single-particle nearest-neighbor hopping. Staying at half-band filling and a 
constant value of U, it was previously shown with Quantum-Monte-Carlo (QMC) 
simulations that one can generate a quantum transition as a function of the coupling 
strength, W, between an AF Mott insulating state and a d x 2_ y 2 SC state. Here 
we complement these earlier QMC simulations with physically more transparent 
diagrammatic calculations. We start with a standard Hartree-Fock (HF) calculation 
to capture the "high-energy" physics of the t — U — W model. Then, we derive and 
solve the Bethe-Salpeter equation, i. e. we account for fluctuation effects within the 
time-dependent HF or generalized RPA scheme. Spin- and charge-susceptibility as 
well as the effective interaction vertex are calculated and systematically compared 
with QMC data. Finally, the corresponding BCS gap equation obtained for this 
effective interaction is solved. 



I. INTRODUCTION 

One salient aspect of the high-T c materials is the vicinity of two, at first sight rather different, 
states of matter, superconductivity (SC) and antiferromagnetism (AF) in their phase diagram. 
The transition between the undoped (AF) system at half-filling and the SC phase is driven by 
doping with mobile holes. In most of the materials, this transition is not direct, and a disordered 
"spin-glass" phase occurs in between. However, it has been argued that the "clean" material would 
display a direct transition from AF to SC phases, and that the spin-glass phase occurs due to the 
high sensitivity to impurity disorder in the vicinity of the phase transition. 

A direct transition from an insulating into a SC phase in a quasi-two-dimensional system (such 
as the high-T c materials) is a very interesting, yet insufficiently understood issue. In fact, it is not 
clear whether this transition is second order down to zero temperature, and thus is related to a 
quantum critical point, or whether there is a finite-temperature classical bicritical point. In the 
framework of the projected SO(5)-theory of high-T c superconductivity!!], it has been suggested that 
the AF and the SC phases may indeed coexist in some portion of the temperature versus doping 
phase diagram. Another open question is the nature, i. e. the universality class of this transition. 
For example, it has been suggested that this transition may be controlled hv an SO(5)-symmetric 
fixed pointcl. SO(5) symmetry is thus restored in the long-wavelength limittl, and AF and SC can 
be described in terms of a unique superspin vectora in the vicinity of the critical point. 

Many efforts have been directed towards studying the AF-SC transition in strongly correlated 
lattice models by numerical techniques such as Quantum-Monte-Carlo (QMC) simulations. As 
a relevant model, the Hubbard model, is widely accepted for the description of salient features 
of high-T c materials. Unfortunately, it is quite difficult to study large enough Hubbard-model 
systems by QMC, due to the occurrence of the minus-sign problem at finite doping. The numerical 
problem can, in principle, be cured, if one can drive the AF-SC transition by means of a parameter, 
alternative to the doping, which conserves particle-hole symmetry and, therefore, avoids the tedious 
minus-sign problem. This idea was followed through by Assaad, Imada, and Scalapino (AIS)EI in 
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terms of their so-called t — U—W model. It rests on adding an interaction term W, which depends 
on the square of the nearest-neighbor hopping. This W term can be obtained from a Su-Schrie 
Heeger type of electron-phonon interaction in the antiadiabatic limitfl. In QMC simulations 
this t—U — W model exhibits a transition from an antiferromagnet to a d-wave superconductor 
at half-filling and at a critical value of the interaction W c w 0.3i (U = At, T = OA). The QMC 
data of AIS supports the picture of a continuous quantum phase transition in the sense that the 
magnetization vanishes continuously at the critical point. The disadvantage of the t — U—W model 
is that the bandwidth grows substantially with W. Therefore, one of usEJ suggested to introduce a 
phase factor in the VF-term which has a d-wave like symmetry. Although this latter model solves 
the problem of the bandwidth, the existence of a phase transition to a d-wave superconductor 
remains open. 

While QMC calculations provide an essentially exact description of the properties of the model, 
semi-analytical, i. e. diagrammatic, calculations allow for a more direct understanding of the pro- 
cesses which are responsible for a given phenomenon. For this reason in this paper, we carry out 
a systematic diagrammatic study of the t—U — W model. 



We first consider in Sec. [II the simple Hartree-Fock level, which, due to the complexity of the 



interaction terms, allows for different broken symmetry phases. However, a careful comparison of 
the energies of these phases shows that the antiferromagnetic phase is always the stablest one, even 
for very large values of W. This holds for both versions of the t—U—W model which are considered, 
i. e. with and without phase factors. Moreover, the only allowed superconducting solution in the 
simple t — U — W model has an s-wave symmetry, while d-wave symmetry is not allowed. These 
mean-field results are in strong contrast with the QMC calculations, which predict a transition to 
a d-wave SC state at some finite WB. On the other hand, in the AF region our mean-field results 
are in very good accord with QMC, in particular concerning single-particle dispersions, as shown 
in Sec. glf 

The fact that the transition to the superconducting state does not come out correctly is of 
no surprise within an Hartree-Fock approximation. Indeed the relevant transition to the d-wave 
SC state is dominantly driven by an effective attractive interaction mediated by spin fluctuation, 
as is the case for the Hubbard modetJo For this reason, one has to consider the effect of 
spin fluctuations beyond the Hartree-Fock level, in order to reach the SC state. This is done 
in Sec. IV, where we carry out a complete RPA summation of all particle- hole diagrams (both 
bubbles and ladders), with Hartree-Fock Green's functions, in order to obtain the frequency- and 
momentum-dependent spin and charge susceptibility. The solution of the corresponding Bethe- 
Salpeter equation is technically quite difficult to. ach ieve and significantly more demanding than 
the standard case of the simple Hubbard modeHEj. This is due to the finite extension of the 
interaction, as well as its dependence on all (three) momenta, and not on the momentum transfer 
only. By changing to a mixed real-space momentum-space representation, we demonstrate that it 
can be reduced, for generic momenta, to the inversion of a 52 x 52 matrix. 

Next, inj-aome analogy to the RPA-analysis of the t — U Hubbard model by Schrieffer, Wen, 
and ZhangE3, we derive the effective two-particle interaction vertex in the static limit, and solve 
the associated BCS equation. As for the simple Hubbard model studied in Ref. [jj], the d-wave 
solution turns out to be the only stable one, in accordance with QMC results. On the other hand, 
as expected, the s-wave solution is unstable, due to the strong on-site Hubbard repulsion U. In 
the d-wave phase, we obtain a decreasing superconducting gap as a function of W, in spite of the 
fact that the attraction between the quasiparticles should be increased by W. This is due, on the 
one hand, to the approximation of taking an energy cutoff for the effective interaction, which has 
been chosen to be of the order of the AF gap, which, in turn, decreases with increasing W. On 
the other hand, the reduction of the density of states at the Fermi level, which is related to the 
broadening of the bands produced by W, contributes in reducing the superconducting gap. 

Our paper is organized as follows. In Sec. || the t — U — W model with and without phase 



factors is introduced, and briefly summarized In Sec. HI, we carry out the Hartree-Fock mean- 
field study of the antiferromagnetic phase. We discuss the HF results and compare them with 
QMC calculations. In Sec. jv] we derive and solve the Bethe-Salpeter equation, i. e. we account 
for fluctuation effects in the time-dependent HF or generalized RPA scheme. We obtain spin and 
charge susceptibilities, as well as the effective interaction vertex. In Sec. [v|, we write down and 
solve the BCS gap equation, obtained from this effective interaction within a static approximation. 



Finally, we present our conclusions in Sec. VI, partly based on detailed comparisons with QMC 
data. 
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II. MODEL 



The Hamiltonian of the t — U — W model is given b> 
with the hopping kinetic energy 

K x = V(ct c- - +4 ? c ) (2) 

<t,<5 

and 



K z = V/(£)(ct c, ? +ct - c ), (3) 
a,Z 

where in its original formulation f(5) = 2. As mentioned in the Introduction, this model was intro- 
duced by Assaad, Imada, and ScalapinotJ, in order torSfeiidy-ithe antiferromagnetic-superconducting 
transition at half-filling. As stressed by these authorsal3'Bll3, the particular choice of W was mainly 
motivated formally as a means of introducing the desired quantum transition from the insulating 
(AF) to the SC state. The choice of the interaction also guaranteed that no fermion-sign problem 
was encountered in the QMC simulations at half-filling. While there are various approximate ways 
to physically justify the form of the microscopic Hamiltonian, it runs into problems when directly 
compared to the high-T c cuprates. One of the problems of this model is the presence of unphysi- 
cally broad bands, which also cause problems in the numerical QMC evaluation. For this reason, 
one of usEJ suggested to introduce in the i^-term of Eq. (0) a d-wave like phase factor of the form 



f(s) = r: t ~ * ■ (4) 




The phase factor restores the correct width of the single-particle bandS|-at half-filling. Unfortu- 
nately, however no transition to a superconductor in QMC simulations!]! has been observed so 
far. 

The W term contains four different processes (see e.g. Ref. [7]), among them single-particle terms 
that renormalize the chemical potential and permit single-particle hopping between second and 
third nearest-neighbor sites as well as singlet and .-triplet scattering terms. These terms are not 
expected to be relevant for the low-energy physical However, the most interesting term for the 
quantum phase transition is the fourth term which generates singlet pair-hopping and produces an 
antiferromagnetic exchange interaction, i. e. 

H® = -2W ]T /(5)/(5')A^,A y , (5) 

7,5,5' 

where At - = (d^d - — cl ct - )/a/2- For 5 = 5' the terms in T-&S contribute to the exchange 

i,5 v i,T i+5,l i,l i+5,l" v w ° 

giving: 

2WY,{S r S u ^\n X n us ). (6) 

i,5 



III. HARTREE-FOCK CALCULATIONS 



The details of our HF calculation are given in Appendix [a|. A fter s olvi ng th e self-consistent 
equations for the mean- field parameters in Eqs. (A5) to ( JA^ ) and (A2C) to (A22) (Appendix [A]) , 
we arrive at the following results: 
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Figure |l| displays the free energy of the different phases (antiferromagnetic, superconducting and 
paramagnetic) for the simple t — U—W model (top) and the t—U — W model with phase factors 
(bottom) as a function of W for fixed U = At and T — OK. For the sake of comparison, we only 
plot the difference to the paramagnetic energy. 

In the simple t — U — W model (Fig. |l|, top) the antiferromagnetic solution (AF) is always the 
most favorable. However, with increasing W the energy of this solution approaches the param- 
agnetic solution (PM). The superconducting solution (SC) has a much higher energy than the 
other solutions. The only possible superconducting solutions have Si = S3 = (see Appendix [a|: 
no on-site pairing, due to U) and s 2 7^ (nearest-neighbor singlet pairing), while for W < 0.3t 
there exists no superconducting solution. The superconducting order parameter s 2 corresponds 
to an s-wave like symmetry. The transition from an antiferromagnet to a (^.^-superconductor 
observed in QMC simulations is not reproduced at the mean-field level. 

In the t — U — W model with phase factors (Fig. [I], bottom) the mean-field ground state is also 
antiferromagnetic (AF). Here, however, in contrast to the simple t—U — W model, the difference 
in energy with the paramagnetic (PM) solution is increasing with increasing W. As discussed 
in Appendix A, there also exist two different superconducting solutions, that lie energetically 
between the antiferromagnetic and the paramagnetic solution and evolve continuously from the 
paramagnetic solution at W = Ot. In the first solution, s\ and S3 are nonvanishing, while S2 = 
(s-wave). The order parameter has a s-wave symmetry with a superimposed weak modulation of 
the gap. In the second, energetically more favorable, solution one has si = S3 = 0, and s 2 7^ 0, 
yielding an order parameter with d-wave symmetry. Notice that also in QMC simulations at half- 
filling, no transition to a superconductor was found in the t — U — W model with phase factors, in 
agreement with our results. However, the mean-field result in Fig. [l] is promising in direction of 
doping away from half-filling, where the AF phase is suppressed. 

The band structure of the antiferromagnetic solutions is evaluated along the usual paths through 
the Brillouin zone, as shown in Fig. ||. Fig. ^ (top) gives the bands of the simple t — U — W model 
for W = 0.15i. One can recognize easily that the bands are much wider than in the Hubbard 
model but their shape is nearly unaltered. This means that if one would scale the bands by a 
factor i, they would be almost identical. The effect of W seems thus to be a mere "dilatation" of 
the bands. 

In the t — U — W model with phase factors, things look quite different. The bands are plotted 
in Fig. (middle) along the path shown in Fig. |](a) and in Fig. || (bottom) along the path shown 
in Fig. Wb) with W = 0.05i. The width of the bands is nearly the same as in the Hubbard 
model, except for the lifting of the degeneracy along the boundaries of the magnetic Brillouin zone 
(MBZ). At k = (it, 0) a kind of double- hump structure can be seen like it appears in t-t'-t" models 
to describe high-T c superconductorsEJ. This can be explained jasJollows: The W term contains also 
hopping processes to second and third nearest neighbor siteaJ'EJ which, due to the phase factors, 
have the same sign as in the standard fit parameters t,t',t" which are often used to adjust the 
bands to the experimental data of Bi 2 Sr 2 CaCu208+5 and YBa2CtePv_s. 

If one compares the antiferromagnetic bands to the QMC datali3'tlfE3 of the t — U — W model 
as in Figs. I and |, one gets a very good agreement. The width of the bands as well as the 
antiferromagnetic gap were reproduced excellently. In addition, the energy bands of the t — U — W 
model with phase factors (Fig. |3|) show the same double-hump structure at k = (it, 0) like it is seen 
in the QMC spectral weight A{k,ui) (Fig. ||). They also reproduce well the lift of the degeneracy 
along the boundaries of the MBZ. 

Finally, we want to look at two characteristic features of the antiferromagnetic solution: the 
sublattice magnetization m and the Mott-Hubbard gap. The sublattice magnetization defined as 

m = \(4^)-(4 d%l )\, (7) 

is plotted in Fig. ^| (top) as a function of W. As expected from the behavior of the free energy 
(Fig. |l|), the sublattice magnetization decreases with increasing W in the simple t — U — W model. 
On the other hand, the sublattice magnetization of the t—U—W model with phase factors is getting 
stronger with increasing W. This is also confirmed by QMC dataHl, which show an amplification 
of the antiferromagnetic correlations with increasing W. 

A similar picture occurs for the antiferromagnetic gap (Fig. g, bottom). Like the sublattice 
magnetization, the Mott-Hubbard gap is decreasing with increasing W in the simple t—U — W 
model, while it increases (nearly linear) with W in the model with phase factors. 
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In summary, the Hartree-Fock calculation gives the antiferromagnetic solution as ground state for 
any values of W in both models. However, qualitatively there are remarkable differences between 
the two models with and without phase factors. A comparison of these results with the QMC data 
shows that the antiferromagnetic phase is described quite well by the mean-field approximation. 

On the other hand, the mean-field level is not able to reproduce the transition to a d x 2_ y 2- 
superconductor at W c » 0.3t observed in QMC simulations in the simple t—U—W model. This 
is of no surprise since results of an Hartree-Fock approximation at finite values of the interaction 
should be taken with due care and cannot give decisive conclusions about the correct phase diagram 
of a model without a comparison with more reliable calculations, such as, e. g. QMC. 

Nevertheless, for large W, for which the AF gap becomes small, one would expect the antifer- 
romagnetic solution to become instable with respect to fluctuations beyond the mean-field level. 
This is what we analyze in the next section. 



IV. TIME-DEPENDENT HARTREE-FOCK (GENERALIZED RPA) 

As demonstrated in the previous section, the HF mean-field approximation is not sufficient 
to describe the transition to a d :E 2_ y 2-superconductor occurring in the simple t — U — W model 
according to QMC simulations. For that reason, we carried out an improved calculation, including 
charge- and spin-density fluctuations. This has been done by means of a time-dependent HF or 
generalized random phase approximation (RPA), in which wc summed both "bubble" and "ladder" 
particle-hole diagrams. In contrast to the fluctuation exchange approximation (FLEX), the Green's 
functions are not calculated selfconsistently, but taken over from the Hartree-Fock results, as it 
has been done in Ref. E3. 



A. Hartree-Fock Correlation Function L° and Interaction Vertex r° 



The 2x2 antiferromagnetic Hartree-Fock Green's function can be written as (see Appendix [A|, 
Eqs. (IaJ) to (|Al|)): 



/ ito + e(k) erA(fc) 

crA(fc) iuj - e(k)J _u,2 _ E 2 (k) ' 

With this Green's function we can construct the Hartree-Fock two-particle propagator L°: 



(8) 



L°o-7C 2 (fcl, fa, Q,Vl,UJ2,v) = 5a 1 <7 2 5a> 1 fc, <Li lo 2 Q H " m2 mi (fa , ^1 , <*\ )Q H " m' x rn> 2 (fa ~ <f, U)\ - V,^), 

"l "2 

(9) 

where the set of m, stand for the indices of the 2x2 matrix in Eq. (||). After a unitary transfor- 
mation with help of the Pauli matrices, i. e. 



where 



L" = U^L°U, 



(10) 

(ii) 



we can write the correlation function in the charge-/ spin-channel representation as: 
Kb™' 1 m ' 2 (fa, fc 2, q) = 



, mi mj 

' L° 00 m i m 2(k u k 2 ,q) 



_ mi rn-2 

\L° z0 «(fa,k2,q) 



<(ki,k 2 ,q) 



L'L + ^<(fa,fa,q) 



L° nz m i m 2(k 1 ,k 2 ,q) 



L° Z z m ' im ' 2 (fa,fa,q)) 

(12) 
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Here and in the following: k = (k,oj), q = (q, v) and so on. 

In contrast to the Hubbard model, also the non-diagonal elements which couple the charge chan- 
nel to the longitudinal spin channel have to be taken into account. For the following calculations 
it is also advantageous to transform from the representation 



to the representation 

L° ab (k 1 ,k 2 ,uj 1 ,uj 2 ;q + nQ,q + n'Q,iy) with 



k u k 2 ,q£MBZ (13) 
ki,k 2 eBZ, qeMBZ, (14) 



where n, n! take the values {0, 1} and Q — 

In this representation, e. g., the longitudinal spin correlation function can be written as a matrix 
in the indices n, n': 



L zz (ki,k 2 ;q + nQ,q + n'Q,v) = 

^ife^ii ( fc i)^n + k 2 +Q§ji ( k i)S 12 i k i -q) \ 

v ^jf^ll (*l)^22(*!l - l) + S k 1 k 2 +Q^ 2 l( k l)G " li k l ~ l) J 

(15) 



Here the Q are spin independent Green's functions (i. e. Q HF (Eq. (@)) with spin index a set 
equal to +1). The interaction vertex in this representation is given by 



Kb( k i, k 2;q + nQ,q + n'Q) = —{ 



1 
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I o) V w {k 2 -q,k u q) 



I ^ 1 )V w (k 2 -q,k\,q + Q)} 



I \)V w {k 2 -q,k\,k\-k 2 )}. 



(16) 



It is written as a direct product of spin- and (n, n')-matrices with the VF-dependent interaction 
V w given by 



V w (k, k', q) = — 8W[(cos fc^, ± cosfcy)(cos k x ± cos k y ) 

+(cos(k' x — q x ) ± cos(fc^ — q y )){cos k x ± cos k y ) 
+(cos(fc s + (fc) ± cos(k y + q y ))(cos k' x ± cosfc^) 
+(cos(fc x + q x ) ± COs{k v + q y ))(cos(k' x - q x ) ± cos(fc^ - q y ))}. 



(17) 



B. Bethe-Salpeter Equation 

With help of the HF two-particle propagator L and the interaction vertex T, we can now write 
down the Bethe-Salpeter equation in the form: 

L ab (k 1 ,k 2 ;q,q',v) =L° a b (ki, k 2 ;£[,{?, v) 

+ ^2 Xic( k i' k 3iq^',v) r ld{ k 3, k 4^\tf")L db (k 4l k 2 ;tf'\c?,v). 

q" ,q"' 
c,d 

(18) 



G 



Unlike for the standard (W = 0) Hubbard model, this equation cannot be easily inverted due to 
the complicated space and spin structure of the W term. For this reason, we apply a method due 
to Hanke and ShamO, which is based on the partial transformation of the interaction vertex back 
into real space. For short-range interactions, this yields finite (small-sized) matrices in real space. 
With the Fourier-transformed correlation function and interaction vertex 

L ah {R x ,R 2 -q,if,v) = -^ e ik ^e- i ^L ab (k 1 ,k 2 ;q^,u), (19) 

ki ,Ai2 

f° ab (R 1 ,R 2 ;q,if) =-JfJ2 e^e-^Tl.ihM-q,^), (20) 

fei ,k 2 

we obtain the Bethe-Salpeter equation in matrix form: 

L ab (Ri,R 2 ;q,cf, v) =L ab (R x ,R 2 ;q,c{ ' ,v) 

+ L° c (i?i, i? 3 ; q, q", v)?° cd (R 3 , R 4 ; q", q"')L db (R A , R 2 ; q"\ if, v). 

This is now a simple matrix equation which can easily be inverted for the interacting two-particle 
propagator L: 

L = (l - L°f "Y 1 L°. (22) 

In contrast to the Hubbard model we have to deal with complex 26 x 26 matrices for the transverse 
spin-channel and with complex 52 x 52 matrices for the coupled charge-/ longitudinal spin-channel. 
The RPA susceptibilities can be constructed by taking the (0, 0)-matrix element in real space, i. e.: 

Xab{q,4\v) = L ab (0,0; q,(f, v) . (23) 

From this one obtains the retarded susceptibilities by the analytic continuation iv — ► lo + ir\. 

Starting from the idea that spin fluctuations are responsible for the pairing of the quasi-particles, 
it is reasonable to first concentrate on the dynamic spin susceptibilities for the antiferromagnetic 
nesting vector Q = (it, n) as a function of oj, as we expect the strongest response there. In Fig. |7] 
(top), Xzz is plotted for the simple t — U—W model with W = O.lt, while in Fig. [t] (middle) Xzz 
is displayed for the model with phase factors and W = 0.05t. Both calculations can be compared 
with the result for the Hubbard model, reported in Fig. [t] (bottom). 

One can clearly see that the spectral weight is mainly concentrated at low frequencies and that 
it is abruptly decreasing at a frequency oj ~ 2A AF , which corresponds to the antiferromagnetic 
gap. This behavior is most evident in the simple t — U — W model. Moreover, one can recognize 
that the overall magnitude of the longitudinal spin susceptibility is biggest in the simple t—U—W 
model and smallest in the t—U—W model with phase factors. 

C. Effective Interaction 

In this section, we calculate the effective 2-particle interaction mediated by the collective charge 
and spin fluctuations evaluated in the preceding section. Here, we restrict to the model without 
phase factor, since this is the only one which, according to QMC calculations, displays <i-wave 
superconductivity. As a first step, we evaluate the fluctuation vertex, from which we can determine 
the modifications of the bare 2-particle interaction given by the t—U—W Hamiltonian. The 
calculation of the fluctuation vertex is performed with the same techniques that were used to 
calculate the Bethe-Salpeter equation. This gives the expression: 
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F a b{ki,k 2 ;q,(f,v) = ^2 T lc(ki,k3;q,tf')L cd (k 3 ,k 4 ;c?\cf'\v)r" db (k 4 ,k2;c?'\cf) 

k 3 ,k 4 

q ,q 

c,d 



E E (w) e ~ iilAi r°ac(Ru R3; q, f) ( 24 ) 



L cd (i? 3 , i? 4 ; 9", 9"', ^)r° db (i? 4 , i? 2 ; <f , ?)e 



which is diagrammatically represented in Fig. |^. 

Next, we have to change from the charge/spin channel representation back to the simple spin 
representation by inverting the transformation given by Eq. (|l0|) : 

+ T z0 (k 1 ,k 2 ;q,q, v)a z ai ff , ct°, a2 +T Z z (k!,k 2 ; q,q , v)o z ai cr*, (T2 

+ r + _(fe,fe;f,^,i/)or- a ,or+ ^ +r_ + (fc 1) fc 2 ;g,g',^+ (7 ,5 : -, (72 ). 

(25) 

Since we want to use the effective interaction in order to write an effective Hamiltonian, Qsly the 
static limit of the fluctuation vertex has to be considered. Thus, simple diagrammatic rulesta yield 
for the correction to the bare interaction: 

(71(72 (71(72 

V °i {h,k 2 ;q, q')=r°i ^(k u k 2 ;q, <f , Q) ■ (-/37V) ■ (-1 ). (26) 
The effective interaction can then be written as (see Fig. ^ for diagrammatic representation) 

V eff °>2(k l7 k 2 ;q,q') = V °i "4 *a; # + V CT ^ ^ (*i , fc 2 ; ?, <f) (27) 



with the bare 2-particle interaction 

(71 (72 

F CT X(fci,fc 2 ;£) =C/<5 CTl(7 ^ CT24 <y ( 7 1 a 3 + F M '(fe2-g,fc 1 ;g)'J < 7 1< 7'/ < 7 2 ^ (28) 

and V w given by Eq. (|l7j). Since one has to consider the pairing of the Hartree-Fock quasiparticles, 
the effective interaction has to be transformed into the 7-base, which produces additional coherence 
factoriM 

For physical reasons, only the pairing of particles with opposite spin (singlet-pairing) was consid- 
ered. In order to take into account the effect of the dynamics on top of our static approximation, 
we follow Ref. [l2|, and introduce a cutoff frequency cj c , analogous to the Debye frequency uj d in 
the standard BCS theory. This ensures that only particles within an interval of width hw c above 
and below the Fermi energy E F are paired. 

The motivation for this cutoff frequency lo c becomes clear if one looks at the spin susceptibilities 
of the simple t—U — W model in Fig. |7[ We have already shown in the preceding section that the 
spectral weight is concentrated at low frequencies. Under the condition that the spin fluctuations 
are responsible for the pairing of the quasi-particles, the longitudinal spin susceptibility gives quite 
naturally a cutoff frequency of the size of the antiferromagnetic gap (uj c ss 2A AF ). 

This implies that for hole dopings away from half-filling, only the intra-valcnce band matrix 
elements have to be considered. The pairing part of the effective Hamiltonian can thus formally 
be written as 

HP " r = E' fc')0(^c - \E°$) - E P \)®{u c - \E°$) - EA)l v l a H% _ a H\_ a ll a , 

k,k' 

(T,(T 

(29) 
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where E v (k) = —E(k) is the valance band energy. 

The direct interaction V^ T {k,k'), which is given by a = a' spin indices contains, besides the 
longitudinal spin fluctuations, also the bare interactions and the charge fluctuations. The exchange 
interaction with a = — a' consists of the transverse spin fluctuations only. In Figs. [To] to O, the 
direct interaction and the exchange interaction were plotted for different paths of k and k' through 
the magnetic Brillouin zone (MBZ) for the simple t — U — W model and the Hubbard model, 
respectively. The exchange interaction was plotted there as — V T"p (fc, —k'), since it has exactly 
this form, with negative sign, in the BCS gap equation (see Eq. (|3l|)). 

Comparing the graphs for the simple t—U—W model (W — O.lt) with the simple Hubbard model, 
one can easily see that the W term amplifies the attractive parts of the direct interaction, whereas 
the attractive parts of the exchange interaction remain constant (see e. g. Fig. |ll|). However, 
the repulsive parts of the direct interaction and the exchange interaction were both attenuated 
considerably by increasing W (see e. g. Fig. |l2|). Therefore, the pairing of the quasi-particles is 
favored in the simple t—U—W model altogether. 



V. BCS GAP EQUATION 

Finally, we want to solve the BCS gap equation for the effective pairing interaction. Starting 
point is the effective Hamiltonian, as obtained in the previous Section 

«e//=£V(£)-M)7t7t 

k . a 

+ 2^ £' V Z%& %W U ° - I 25 "® - E r \)Q( Uc - \E»{k>) - EAhl)^% t _^U,-A- 
%,%' 

I 

(30) 

With this Hamiltonian we want to study the superconducting properties of the simple t—U — W 
model for different hole doping and different values of the model parameter W . The BCS gap 
equation becomes 

A & = "4 EVn& *') - V n (k, -k'))^X. (31) 

Here we used the abbreviations 

E{k) = y/?{k) + A*(k), (32) 

(33) 

V acr ,{k,k') = V^j{k,k')Q{u A - \E v {k)-E F \)Q{u A - \E"{k') - E F \). (34) 

The gap equation ( ]3l] ) was iterated by assuming different symmetries of the superconducting order 
parameter, however, only d-wave solutions turn out to converge. 

The results for the simple t—U—W model are shown in Fig. |l3| as diamonds, triangles and 
squares. For all values of W considered, the superconducting gap becomes zero at half filling, 
which is consistent with QMC results. It can be shown easily with the aid of equation ([H]), that 
even including interband matrix elements doesn't change this result. This is due to the fact that 
the band gap is still quite large for these values of W (cf. Fig. ||) for interband matrix elements 
to contribute substantially to the gap equation. On the other hand, all curves seem to indicate 
that the superconducting phase starts at very small doping, which is in contrast with experiments. 
While there are no conclusive QMC results about the t — U — W model in this region, due to the 
minus sign problem, the limitations of our perturbative procedure applied for moderate values of U 
suggest to consider this result with due care. Indeed, strong phase fluctuations, not included in a 
BCS-type of calculation like Eq.(|3~l]) are known to be important and to suppress superconductivity 
at small doping. 
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Figure |13| also shows, that the optimal hole doping for the simple t — U — W model is moving 
closer to half-filling with increasing W. For comparison, we show in Fig. |l^ (stars) a result for 
the t — U — W model with phase factor. One can see, that the superconducting gap is strongly 
suppressed for hole doping near half-filling, compared to the simple t — U — W (and Hubbard) 
model. The effect is here even stronger than for the simple t — U — W model, since in the model 
with phase factors the gap increases with W. 

It remains the question as to why the superconducting gap is getting smaller with increasing W in 
the simple t — U — W model, in spite of the fact that the bare attraction between the quasi-particles 
is enhanced by the W term. To understand this, one has to keep in mind that there are other 
important quantities, like the cutoff frequency or the density of states at the Fermi level, which 
have an important effect on the magnitude of the superconducting gap. How the superconducting 
gap depends on these quantities is qualitatively seen already in the weak-coupling solution of the 
original BCS equation for an attractive 5-potentialtLj: 

A = 2^ D cxp(-^-). (35) 

A stronger coupling g is increasing the superconducting gap, while a smaller cutoff frequency 
uj d is decreasing it. Moreover, if the density of states at the Fermi level iV(0) is reduced, the 
superconducting gap is getting smaller, as well. This final point is decisive, since the W term 
is broadening the energy bands considerably, thus reducing the density of states dramatically. 
Therefore, the fact that the superconducting gap is getting smaller with increasing W can be 
understood within our approximation. 



VI. COMPARISON WITH QMC RESULTS AND CONCLUSIONS 

In this paper, we have first shown that the standard Hartree-Fock approximation can describe 
the antiferromagnetic properties of the t — U — W model, especially the single-particle energy 
bands, in surprisingly good agreement with the QMC simulations. On the other hand, like for the 
Hubbard modelE3, it is unable to reproduce the transition to a d-wave superconductor observed in 
the simple t — U — W model in numerical simulations (QMC). In order to partly overcome this 
short-coming we adopt a time-dependent HF or generalized RPA calculations, which we present 
in the second part of our paper. 

The standard Hartree-Fock approximation captures only the "high-energy" physics, and is thus 
capable of reproducing band- widths and the overall features of the single-particle spectral-function. 
Equivalently, short-range pairing-correlation functions should be well reproduced within this ap- 
proximation. This is indeed the case. At short length scales (as shown in table 1) the extended 
s-wave vertex contribution to the pairing correlation functions is dominant in QMC simulations of 
the simple t-U-W model. It is only at larger distances that the d-wave pairing correlations become 
dominant. This crossover from short range to long range properties is not reproduced within the 
standard mean-field approximation. Alternatively, for the t-U-W model with phase factors d-wave 
pairing dominates in QMC simulations at small length scales (see table 2). 

On the other hand, by including charge and spin fluctuations within a time-dependent HF, or 
generalized RPA summation of ladder and bubble diagrams, we were able to obtain an effective at- 
traction between the quasi-particles, which is enhanced by W. Moreover we obtain a corresponding 
superconducting order parameter with the correct d-wave symmetry. 

The QMC results indeed show that as W increases at fixed U or U decreases at fixed W, an 
instability towards d-wave superconductivity occurs. To illustrate this, Fig. |lj plots the vertex 
contribution to the <i-wave pairing correlations as well as the staggered spin susceptibility. As ap- 
parent at low U for fixed W, the superconducting d-wave becomes the leading instability. Finally, 
it should be pointed out that our Hartree-Fock-Bethe-Salpeter procedure is of perturbative nature 
and, therefore, our results should be considered with due care, due to the fact that the interaction 
is not small. 
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APPENDIX A: DETAILS OF THE HF CALCULATIONS 



1. Antiferromagnetic Mean Field 

We start with the following Ansatz for the mean field parameters: 

( C L C ?J = n 1 +ae i & i n 2 , 



i+S',tr' 



n 5 , 



(Al) 
(A2) 
(A3) 



where Q denotes the antiferromagnetic nesting- vector Q = (ir,n). Here the expectation values 
( • • • ) contain also an average over all 8 and 5' , whenever explicitly present. At half filling, it is 
easy to show that n\ = ^, and 114 = |. The mean-field Hamiltonian can then be written as 



nj M F = Sr'(^ ?t N ( £ (k) <jA(k)\ I 



+ Et 



with 



and 



e{k) = — 2i(cos k x + cos k y ) — 96VFn 3 (cos k x ± cos k y ), 
A(k) = -Un 2 + 32Wn 5 - 8Wn 2 (cosk x ± cosfc^) 2 

Kuw = +UNnj + 192WNnl - 64WN{± + n 2 n 5 ). 



(A4) 

(A5) 
(A6) 

(A7) 



Here, and in the following equations, the upper sign stands for the simple t—U — W model, while 
the lower sign is used for the t—U — W model with phase factors. The Hamiltonian can be 
diagonalized by the usual transformation 



u{k) <jv{k) 
v(k) —au(k) 



with 



u(k) 

v(k) 



1 + 



e{k) 
E(k) / 

E(k), 



and 



E(k) = y / e 2 (fc) + A 2 (fc). 



(A8) 

(A9) 
(A10) 

(All) 
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The resulting Hamiltonian is given by 



h mf = X^' E(k)h? 7 £ -7^ 7 H ) + E M 

tl -tuw V ^ 'fc,<T 'fe,<T 'fc.cr ik,<r> T tCJ 



MF 
W ' 



(A12) 



k.a 



Already at this point, one can see from equations (A5) and (All) that the parameter n 3 produces 
extremely wide bands in the simple t — U — W model. On the other hand, it is straightforward to 
see that for the alternative t—U — W model with phase factors one must have n 3 = 0, in order to 
preserve the symmetry of the energy bands under interchange of x- and y-directions. This explains 
why the bandwidth of the t — U — W model with phase factors is drastically smaller, and essentially 
the same as the one of the simple Hubbard model. 



2. Superconducting Mean Field 

Here, the mean field parameters are chosen as: 



(ct c, ) 


= ni, 


(A13) 


(ct f{6)c~ x ) 


= n 3 , 


(A14) 




= "-4, 


(A15) 


(ct ct ) 

i,a i, — cr 


= CSi, 


(A16) 


( c t f(5)d T ) 


= ers 2 , 


(A17) 


ct ? f{S)f(S')cl , ) 

i+S,a J WJ y ' i+S',-a' 


= cs 3 . 


(A18) 



The superconducting parameters can be divided into two groups: S2 stands for the nearest-neighbor 
singlet pairing, which is favored by the W term (Eq. (|^)). On the other hand, s\ and S3 for 8 = 5' 
represent the on-site singlet pairing which is suppressed by the Hubbard U. The mean-field 
Hamiltonian thus is 

H MF -V(c f c - )( £ $) A ^}] ( %A+E MF (A19) 

n tuw - Z>- T ,c_ M j ^ A( ^ +* tuw (aw) 

with the single-particle energy e(k) and the gap parameter A(fc) given by: 

e(k) = —2t(cosk x + cosk y ) 

-96Wn 3 (coak x ± cos k y ), (A20) 

A(fc) = +Usi 

— 8M^si(cos k x ± cos k y ) 2 

— S2W s 2 {cosk x ± cosfc y ) 

— 32Ws 3 . (A21) 

The energy constant E™ 1 ^, stands for: 

E™w = ~UNs\ + 6AWN( Sl s 3 + s\ - i) + l§2WNn%. (A22) 
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Eq. ( A19) can be diagonalized with a Bogoliubov - de Gennes transformation, similar to Eq. ( A8), 
i. e. 



v(k) -u(k)J IS* 



u{k) 



v{k) 



u{k) v{k) 
%{k 

e{k) 



1 



E{k), 
E{k), 



sign(A(fc)). 



The resulting Hamiltonian now becomes 



MF 
UW > 



(A23) 
(A24) 
(A25) 

(A26) 



with the usual relation (Eq. ( |All )) for E{k). 

The para magnetic solutions can be easily obtained by setti ng the antifcrromagnetic parameters 
in Eq. (A4), or the superconducting parameters in Eq. (A19), equal to zero. 
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APPENDIX: FIGURES AND TABLES 



t-U-W model 



— 

a 




0.1 0.2 0.3 0.4 0.5 
model parameter W/t 

t-U-W model with phase factors 




0.1 0.2 0.3 0.4 0.5 
model parameter W/t 

FIG. 1. Free energy per lattice site of the different phases of the simple t — U — W model (top) and 
the t — U — W model with phase factors (bottom). Plotted is the difference to the paramagnetic energy 
(U = At, T = OK, ii = 1, lattice size: 60 x 60). 
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T * 



(a) (b) 
FIG. 2. Paths followed through the Brillouin zone to plot the energy bands of Fig. |^. 



t-U-W model 




(0,0) 



(71,0) (7C,7C) (71/2,71/2) (0,0) 

wave vector k 



t-U-W model with phase factors 




(0,0) (7C,0) (71,71) (7C/2.7C/2) (0,0) 

wave vector k 

t-U-W model with phase factors 




(it,0) 



(71/2,71/2) 
wave vector k 



(0,71) 



FIG. 3. Energy bands of the simple t — U — W model (upper graph, W = 0.15t) and of the t — U — W 
model with phase factors (middle and lower graph, W = 0.05t) on the paths shown in Figs. |^(a), for upper 
and middle graph, and ^(b) for lower graph (U = it, T = OK, (j, = 0t, lattice size: 100 x 100). 
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FIG. 4. Spectral weight A(k,ui) of the simple t — U — W model, obtained from QMC simulations (U = At, 
W = 0.15i, T = OK, fx = Ot) (from 0). The considered path in the Brillouin zone is listed on the left 
hand side of the figure. We have normalized the raw data by the factor listed on the right hand side of 
the figure. This normalization sets the peak value of A(k,u) to unity for all considered k vectors. 



T = 0, (w) = 1, UJt = 4, W d /t = 0.05, L = 12 




FIG. 5. Spectral weight A(k,uj) of the t — U — W model with phase factors, obtained from QMC 
simulations (U = At, W = 0.05t, T = OK, fi = Ot). 
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o 0.8[ 

N 

u 0.6[ 
a 

5P 



o 



0.4- 



original t-U-W 

t-U-W with phase factors 



| 0.21- 

3 



0.1 0.2 0.3 0.4 0.5 
model parameter W/t 



Mott-Hubbard gap 



1 



original t-U-W 

t-U-W with phase factors 



o 
S 2 



nl i i i i i i i i i i i I 

u 0.1 0.2 0.3 0.4 0.5 
model parameter W/t 

FIG. 6. Sublattice magnetization (top) and Mott-Hubbard gap (bottom) of the t — U — W model as a 
function of W (U = At, T — OK, fi — Ot, lattice size: 100 x 100). 
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1,1,1,1 



12 



co/t 



16 



t-U-W model with phase factors 




co/t 

Hubbard model 




FIG. 7. Longitudinal spin susceptibility x^^iQy Q\ w ) °f the simple t — U — W model (top, W = 0. 
the t — U — W model with phase factors (middle, W = 0.05i) and the Hubbard model (bottom, W = 
(as usual: U = At,T = OK, jj, = 0t). 
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FIG. 8. Diagrammatic representation of the fluctuation vertex 




FIG. 9. Diagrammatic representation of the effective interaction 
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k,k' 

Hubbard model 



> 



v 7(w 



k,k' 

FIG. 10. Direct part V£" r (k, k') and exchange part — V£™^(k, — k') of the pairing-potential for the 
simple t~U — W model (top, W = O.lt) and for the Hubbard model (bottom, W = Ot). The paths shown 
above were followed by k and k' through the MBZ. The distance between k and k' amounts 8k ~ 0.147T 
(U — 4t, T = OK, n, = Ot, lattice size: 100 x 100). 
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FIG. 11. Direct part V*% tr (k,k') and exchange part — V^"^(k, — k') of the pairing-potential for the 
simple t — U — W model (top, W = O.lt) and for the Hubbard model (bottom, W = Ot). The path shown 
above was followed by k' through the MBZ with k held constant (U = 4t, T = OK, fi = Ot, lattice size: 
100 x 100). 
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t-U-W model 




k,k' 

FIG. 12. Direct part V£" r (k,k') and exchange part —V^^(k,—k') of the pairing-potential for the 
simple t — U — W model (top, W = O.li) and for the Hubbard model (bottom, W = Ot) The paths shown 
above were followed by k and k! through the MBZ. The distance between k and k! amounts St « 0.2n 
(U = 4t, T = OK, fx = 0t, lattice size: 100 x 100). 
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FIG. 13. Maximum value of the superconducting gap function A maa . = max A(fc) as a function of the 
mean electron number per lattice site (n) for different values of the model parameter W (U = At, T = OK, 
lattice size: 20 x 20, d-wave symmetry). The curves for the simple t — U — W model are depicted as 
diamonds, triangles and squares and the curve for the t — U — W model with phase factors is depicted as 
stars. 
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FIG. 14. Vertex contribution to the d-wave pairing correlations (•) and staggered spin susceptibility 
(O) for the simple t — U — W model, obtained from QMC simulations. 
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r 


TO 


TO 


(0,0) 


0.2950 ±0.0018 


0.1304 ±0.0011 


(0,1) 


0.0932 ± 0.0009 


0.0238 ± 0.0006 


(0,2) 


0.0076 ± 0.0002 


0.0252 ± 0.0003 



TABLE I. Short range vertex contribution of pair-field correlations in the extended s- and d-wave 
channels for the simple t — U — W model, obtained from QMC simulations. Here we consider an L — 24 
lattice at W/t = 0.35, U/t = 2 and (n) = 1. We choose antiperiodic (periodic) boundary conditions in the 
x (y) direction. The distance r is in units of the lattice constant. 



f 


TO 


TO 


(0,0) 


1.661 ±0.079 


3.082 ± 0.066 


(0,1) 


0.373 ± 0.020 


1.107 ±0.015 


(0,2) 


0.091 ±0.008 


0.119 ±0.004 



TABLE II. Short range vertex contribution of pair-field correlations in the extended s- and d-wave 
channels for the t — U — W model with phase factors, obtained from QMC simulations. Here we consider 
an L — 10 lattice at W/t = 0.35, U/t — 4 and (n) = 1. We choose periodic boundary conditions. The 
distance r is in units of the lattice constant. 
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